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1 Introduction 



The solution of the Schrodinger equation for any spherically symmetric potential has attracted atten- 
tion over the past years [1-16]. The motivation in this direction arises from considerable applications 
in the different fields of the material science and solid physics. 

The Harmonic oscillator and H-atom (Coulombic) problems are two exactly solvable potentials 
which have been thoroughly studied in A^-dimensional space quantum mechanics for any angular mo- 
mentum i. These two problems are related together and hence the resulting second-order differential 
equation has the normalized orthogonal polynomial function solution (cf. Ref.[17] and the references 
therein). On the other hand, the Pseudoharmonic and Mie type potentials are also two exactly solv- 
able potentials other than Coulombic and Harmonic oscillator, their wavefunctions are vanishing at 
the origin. 

The path integral solution for one-dimensional special case Mie-potential which is a perturbed 
Coulombic-type potential was obtained in [18]. Moreover, the Schrodinger equation for a system 
bound by a Mie-type potential was also solved by using the 1 /N expansion method [19]. 

In this letter we will follow parallel solution to Refs.[17, 20] and give a complete normalized 
polynomial solution for the general A^-dimensional Schrodinger equation for diatomic molecular sys- 
tems interacting through Mie type potential, Coulombic-type potential with an additional centrifugal 
potential barrier, which reduces to the standard three dimensional case when the parameter N is set 
equal to 3. Our aim is to present a different approach to calculate the non-zero angular momentum 
solutions of the A^-dimensional space of Schrodinger equation for the Mie-type potential. 

The contents of this paper is as follows. In Section |2l we give the eigensolution of the di- 
mensional Schrodinger equation with Mie type potential for the analytical bound-state (real) and 
imaginary solutions. Finally, in Section[3]we give our results and conclusions. 
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2 The A^-Dimensional Schrodinger Equation with Mie-Type Po- 
tential 



We wish to solve the A^-dimensional Schrodinger equation for a special case of the Mie-type potential 
[18] given by 



V{r) = Do 



k fro 



m ffoy 



(1) 



m — k \ r J m — k \ r J 
where Do is the interaction energy between two atoms in a solid at r = tq, and m > A; is always 
satisfied. By taking m = 2k and choosing the special case k = 1, corresponding to a Coulombic-type 
potential with additional centrifugal potential barrier, Eq.(l) becomes 



N 

2 \ ^ 2 

i=l 



(2) 
(3) 



where oi = 2Doro and 02 = 2Doro^. For brevity, we write the radial part of the Schrodinger equation 
in A/^ -dimensions as 



which reduces into the form 



(4) 



(P N -1 d l{l + N -2) 2m 



I dr'^ r dr 
Using the dimensionless abbreviations: 



Enl H ttt; 



02 



Rni{r) = 0. 



(5) 



X = r/ro;/3^ 



2mrl 2 2mrQ 



gives the following simple form 



(6) 



d^Rnlix) , N-ldRnlix) 



dx'^ 



+ 



X dx 



+ 



2 272 j^ + l(l + N- 2) 



-13' + ^ 

X 



x^ 



Rni{x) = 0, (7) 
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which restores its three-dimensional Schrodinger standard form once N — 3. Eq. (7) has an irregular 
singularity in the x ^ oo limit where its normalizable solutions in bound-states behave as [17, 20] 



{ — -(3^]Rni{x)^0, (8) 



so that Rni{x) — Ani exp(— + Bni exp{Px) In order that Rni{x) — > as x — > oo, we set Bni — 0, 
so 

lim Rni{x) = Aniexp{-Px), (9) 

x— »oo 

where A^i is the normalization constant. Eq. (8) leads us to propose the trial solution 

Rni{x) = Nni ex.p{-/3x)g{x), (10) 

where A^^; is another normalization constant. Putting this back into Eq. (7) yields an equation for 
g{x) of the form 

where the prime refers to the derivative with respect to x. Because exponential behavior has already 
been taken out, one hopes that the solution for g{x) is a polynomial. Indeed, Eq. (7) has a singularity 
at X — > 0, the substitution of the trial solution g{x) — x^, provides the positive root solution: 



2 V V 2 

Asq> 0, the wavefunction vanishes at x = 0, corresponding to the strong repulsion between the two 
atoms. It is reasonable to set 

Rni{x) = N^ixiexp{-Px)h{x), (13) 

to Eq. (7), one gets 

h"{x) + - 2/?) h'{x) 

+ - - (iV - 1)/^ ^ q{q-l) + {N-l)q-^'-l{l + N-2) \ ^ ^ ^^^^ 

\ X x"^ I 



Setting the numerator of term equal to zero, in the last equation, and solving the resulting 
quadratic equation leads again to the solution given in Eq. (12) and thus gives the following dif- 
ferential equation 



xh"{x) + [2q + N-l- 2Px] h'{x) + [27^ - 2g/3 - {N - 1)/?] h{x) = 0. (15) 

The confluent series, for large values of x, is proportional to exp(2/9a;) so that Rni{x) diverges for 
X — > 00 if the series iFi does not break off. If it does iFi is a polynomial and Rni (x) ^ for x — > 00 
becomes normalizable. After the substitution of series form [20] 



h{x) = £ Qx\ (16) 

1=0 



into Eq. (15), it provides 



_ i + g + (iV-l)/2-7V/3 ^ _ „ „ C,+i 1 

- {i + m + 2q + N-l) - ^ ~C~ V ^^^^ 

which leads to a divergent wavefunction if not truncated to a maximum value for i. Neverthless, the 
wavefunction is finite everywhere since i and I are finite, it follows 

iN-l) 72 

Vax + g H ^ ^ ^ ^' ^n;{n = 0, 1, 2, ...) . (18) 

Further, using the following abbreviations in Eq. (15): 

z = 2/3x; c = (g + ; -n = (c/2 - 7V/3) , (19) 

gives the general type of Kummer's (Confluent Hypergeometric) differential equation of the form 

zh"{z) + [c-z] h'{z) + nh{x) = 0; (n = 0, 1, 2, ...) , (20) 

having the solution 

^■^ (22 2n—l) 

h{x) =1 Fi(a, c; z) =1 Fi(-n; 2j- - 2n; 2/3x) = Vn ' \2(3x), (21) 

where iFi(— n, m + 1; 2;) = L^^"^) [z) denotes the Kummer's function. On the other hand, we may also 
write Eq. (15) in the following general form 



zh"{z) + [m + l- z] h'{z) + nh(x) = 0, 



(22) 



where m = 2q + N — 2 and rimax — — Q — ^^2^'^ = 0, 1, 2, .... If n is a non-negative integer, 
then a finite polynomial solution is allowed. (This when combined with the rest of Rni{x) yields a 
normalizable solution.) In particular, this solution to Eq. (22) is the generalized Laguerre polynomial 

L^\2(3x). Combining everything one finally has 



nl 



r \ — 



(N-l) 



exp 



Since the radial volume element in A^-dimensional space is r^~^dr, one can obtain [17] 



nl 



-1/2 



(23) 



N, 



nl 



m 



q+N/2 



(m + n)! 



nl 



{2n + m+ 1), 



or equivalently as 



N, 



nl 



nl 



^r(f -n) 



1/2 



(24) 



(25) 



2.1 Negative Energy 

For bound-states, P > 0, the solution becomes in arbitrary normalization constant is given by means 

of Eq. (23) in which Rni{x) for x — > 00 becomes normalizable. On the other hand, from Eq. (6) 
with the aid of Eqs. (12) and (18), the eigenvalue becomes [21] 



2mrQ 



n + l + J{l + ^)\,^ 



(26) 



Since the parameter 7 S> 1 for most diatomic molecules, we may expand into powers of I/7. This 
leads to 



Do 



7 



+ 



(27) 
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The Mie-type potential, Eq.(2), can be expanded about its minimum atr — as 

v(r) ^ D,^S^ _ n,. 

Therefore, with classical frequency for small harmonic vibrations. 



(28) 



2Dn 



CO — 



2 ' 



(29) 



and the moment of inertia 



mro, 



(30) 



we amve at 



nl 



3^ 
21 



1x2 



3h^ ( 1\ /, N -2\^ 



(31) 



2.2 Positive Energy 



For /3 < or E^i > 0, it is no longer real but purely complex, Px — —inr with n 
gives the wavefunction [21] 



^ which 



N 1 i'y'^ 

Rni{r) = Ar,i{r/roYeMii^r)iF^{q + - - - - 2g + AT - 1; -2mr), (32) 
which vanishes at r = 0. Its asymptotic may be found from the well-known formula [21] 



iFi(a, c; z) ^ exp(-i7ra)-i^z-'^ + ^ exp(z)z"-^' = 0, 

1 c — a la 



(33) 



holding for the whole complex 2;-plane cut along the positive imaginary axis. 



Rni{r) = C„/r''exp(i«;r) 



-f^-^+f-l) r(2g + iV-l) 



-2iKr) 



i7'' (JV-1) 
«rQ 2 



r(2; + Af-i) 



—2iKr) "'■0^ 



(34) 



with^ 

Kro 



9 D"^ 

^-^ro- Therefore, the eigenenergy reads as 



nl 



2mrQ 



n+- + 



(35) 



3 Results and Conclusions 

We have studied the solution for a Mie-type potential. Considering the special case of the Mie poten- 
tial with m = 2k and for k = 1, the problem was reduced to a Coulombic potential with the additional 
centrifugal potential barrier of order The exact solutions for this particular case have been ob- 
tained, which are similar to the Hydrogenic solutions [17]. We have calculated the eigenvalue and the 
corresponding wave function considering the negative (bound-states) and positive (imaginary) cases 
for any quantum-mechanical system bound by a special case of the Mie potential. 

The present results for the potential parameters ai — 1 and 02 = in Eq. (2) correspond to the 
Coulombic potential case. They are the same with the previous calculations in [17-19]. Thus, for this 
particular case the energy terms in the expansion, Eq. (27), take the form 



En — 



(36) 



This value is the exact nonrelativistic H-atom energy expression. Further, setting N = 3, Dq = Vo/2, 
tq — a, A — ^g'^Vq and B — g'^Vq, Eq.(27) becomes 



H/r) — rr-D 



2n' + 1+ (2/ + 1)' + 



ft' 



1/2- 



n = n 



1 = 1,2,3... 



(37) 



which consequently recovers the formula (28) in Ref. [18] and also Eq. (19) in Ref. [19]. For the 

H-atom case this energy expression also gives —0.5 a.u. exactly. Thus, the present result reproduces 
exactly the path integral, 1/A^-expansion, and the non-relativistic Schrodinger equation solutions. 
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